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Abstract 

For the three-dimensional cubic model, the nonlinear susceptibilities of the 
fourth, sixth, and eighth orders are analyzed and the parameters 5^^^ charac- 
terizing their reduced anisotropy are evaluated at the cubic fixed point. In 
the course of this study, the renormalized sextic coupling constants entering 
the small-field equation of state are calculated in the four-loop approxima- 
tion and the universal values of these couplings are estimated by means of 
the Pade-Borel-Leroy resummation of the series obtained. The anisotropy 
parameters are found to be: 6^^'^ = 0.054 ± 0.012, 6^^'^ = 0.102 ± 0.02, and 
^(8) = 0.144 lb 0.04, indicating that the anisotropic (cubic) critical behavior 
predicted by the advanced higher-order renormalization-group analysis should 
be, in principle, visible in physical and computer experiments. 



PACS numbers: 75.40. Cx, 64.60.Ak, ll.lO.Lm 



I. INTRODUCTION 

Three decades ago Wilson and Fisher discovered that in the vicinity of the critical point 
the order parameter fluctuations may change the effective anisotropy of the system drasti- 
cally [jl|. Having studied the anisotropic XY model within the e expansion they found that 
approaching this model either becomes effectively isotropic or developes its anisotropy 
further |l[] until the fluctuation-induced first-order phase transition takes place . What 
scenario the system chooses depends on the value and sign of the anisotropic coupling con- 
stant in the initial Landau- Wilson Hamiltonian. In the course of studying of the generalized 
cubic model, the crucial role of the order parameter dimensionality n was revealed: for 
n < rif. the model undergoing the continuous phase transition demonstrates isotropic critical 
behavior while for n > Uc it remains anisotropic at criticality 0]. The numerical value of 
the marginal spin dimensionality ric separating these two regimes is of prime physical im- 
portance since it fixes the true mode of the critical behavior of real cubic ferromagnets and 
of some other systems of interest. 
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First estimates of Uc for the cubic model were obtained within the lower-order 
renormalizat ion-group (RG) calculations and turned out to be in favor of the inequal- 

ity > 3 indicating that cubic ferromagnets should belong to the class of universality of the 
three-dimensional (3D) Heisenberg model. The classical technique of the high-temperature 
expansions seemed to support this conclusion [Q. However, the resummation of the three- 
loop RG expansions [^[101 and subsequent higher-order RG calculations performed both in 
three ITTl-p^ and (4 — e) (Ref. [Q) dimensions shifted the numerical estimate for down- 
wards, fixing it below 3 [P-|2l[|. Remarkable consensus between different field-theoretical 
approaches and resummation techniques was achieved in the course of this study: 3D RG 
calculations |jl3|, resummed e expansion biased e expansion technique |jl3|, and pseudo-e 
expansion analysis |jl9| lead to ric = 2.89, ric = 2.855, Uc = 2.87, and ric = 2.862, respec- 
tively. Moreover, the predicted numerical value of Uc was found to be stable with respect 
to the order of perturbation theory: starting from the four-loop approximation, the 3D RG 
estimate stays practically unchanged varying around 2.9 within the interval ±0.02 |[TT|-[l3 . 
Hence, according to the up to date, most accurate theoretical data, cubic ferromagnets 
should demonstrate the specific - cubic - critical behavior with a special set of critical 
exponents. 

Since Uc is very close to the physical value n = 3, the fixed point governing the cubic 
critical behavior lies very near the Heisenberg fixed point at the RG flow diagram and the 
critical exponents for both fixed points almost coincide. For instance, the susceptibility 
exponent 7 is equal to 1.3895(50) for the 3D Heisenberg model [^] while at the cubic 
fixed point 7 = 1.390(6) |13|. It is clear that measuring the critical exponents one cannot 
distinguish between the cubic and Heisenberg critical behaviors. 

In such a situation, some alternative physical quantities should be addressed to clear up 
how the systems with the cubic symmetry behave approaching Tc. In this paper, the nonlin- 
ear susceptibilities of the cubic model are studied in the framework of the field-theoretical 
RG approach in three dimensions, with particular attention paid to their anisotropy at crit- 
icality. It will be shown that the anisotropy of nonlinear susceptibilities is sensitive to the 
type of critical asymptotics and its measurement can be used, in principle, for detection 
of the cubic (non-Heisenberg) critical behavior, provided the anisotropy is as strong as the 
modern higher-order RG calculations predict. 



II. NONLINEAR SUSCEPTIBILITIES AND EFFECTIVE COUPLING 

CONSTANTS 

In the critical region, the expansion of the free energy of the cubic model in powers of 
the magnetization components Mq may be written in the form: 

F(M„, m) = F(0, m) + ^m^'W^ + m^'^'^{ui + Vi5^p)MlMl 

+m-3^(M6 + q^5^p + v^5^p5o,^)MlMlM'^ + ... (1) 

where 77 is a Fisher exponent, m being an inverse correlation length, and M4, ^4, Wg, q^-, and 
Vq are dimensionless effective coupling constants acquiring, under T — >• T^, certain universal 
values. These coupling constants are related to the nonlinear susceptibilities defined in a 
conventional way: 
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Of prime importance are the nonlinear susceptibilities in two particular cases, when (i) an 
external magnetic field is parallel to a cubic axis (x^*-*) and (ii) it is oriented along a cell 
space diagonal (x^'*)- For these two highly symmetric directions the differences between the 
corresponding values of the nonlinear susceptibilities are known to be maximal: i. e., the 
anisotropy is most pronounced. It is easy to show that 
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where x is a linear susceptibility. To characterize the anisotropy strength, we define the 
reduced parameters 
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Below, they will be estimated at criticality. 

Let us start from the lower-order nonlinear susceptibility x^^'^ ■ The coordinates of the 
cubic fixed point u*^ and f 4 in three dimensions are known from the higher-order RG calcula- 
tions, with resummed four-, five-, and six-loop RG expansions yielding very close numerical 
Considering recent six-loop RG estimates as the most reliable ones, we accept 
0.755 ±0.010, vl = 0.067 ±0.014 The substitution of these numbers into Eqs. (3) 
and (5) gives 



results ||T1|JI3 



5(4) = 0.054 ±0.012. 



(6) 



In the next section, the universal critical values of the sextic coupling constants will be 
evaluated and an estimate for 6^^^ will be obtained. 



III. RG EXPANSIONS FOR SEXTIC COUPLINGS AND ANISOTROPY OF 

SIXTH-ORDER SUSCEPTIBILITY 

Let us proceed to determination of x''c\ xf\ and 6''^\ As seen from Eq. (4), we are in 
a position to calculate the effective coupling constants uq, qq, and vq. They will be found 
perturbatively, using the field-theoretical RG approach in three dimensions. Our analysis is 
based on the well-known Landau- Wilson Hamiltonian of the 3D n-vector cubic model: 
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where is the reduced deviation from the mean-field transition temperature. All the 
RG calculations are carried out within a massive theory under zero-momentum normaliz- 
ing conditions for the renormalized Green function m) and the four-point vertices 
Unipi, m, u,v,), Vnipi, m, u,v,): 



G^i(0,m) = m^, 
Uji{0,m,u,v) = mu, Vr(0, m, u, f ) = mw. 



gp2 



(8) 



Here, the value of the one-loop vertex graph including the factor {n + 8) is absorbed in u 
and V in order to make the coefficient for the term in the function equal to unity. The 
quartic effective couplings u and v thus defined are related to M4 and entering Eqs. (1), 
(3), and (4) in an obvious way: 
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n + ^ 
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We limit ourselves to the four-loop RG approximation that proved to lead to quite a 
good numerical estimate for the universal value of the sextic effective coupling for the 3D 
0(n)-symmetric model |^,24|. Since the symmetry factors and integrals for all the relevant 



Feynman diagrams have been found earlier ||2^, what we have to do here is the calculation of 
the tensor (field) factors for the six-point vertex graphs and mass insertions generated by the 
0(n)-symmetric and cubic interactions. Performing this calculation and then renormalizing 
the "bare" perturbative expansions for the six-point vertices precisely in the same way as 
has been done for the isotropic n- vector model p3| , |2^ , we arrive at the following four-loop 
RG expansions: 
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+ (19.382741n + 1873.9825)m2^; + (2.3101121n + 679.33934)ui;2 + 100.90652^^ 

- - — ^—^ [(-0.0638n^ + 52.451^2 + 2314.99n + 14387.6)u^ 
+ (-2.3434n2 + 1006.56n + 29601. 0)M^f + (100.156n + 18985.4)mV 
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These RG expansions should obey some exact relations appropriate to systems with 
cubic anisotropy. Indeed, for n = 2 the model (7) is known to possess a specific symmetry 
property. If the field (pa undergoes the transformation 



^1 + ^2 



V2 ' ' 72 

the quartic coupling constants are also transformed, 



u 



3 

U+-V 



(13) 
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but the structure of the Hamiltonian itself remains unchanged [|l| . Since the RG functions of 
the problem are completely determined by the structure of the Hamiltonian, the RG equa- 



tions should be invariant with respect to any transformation conserving this structure p5 



This should be also true for all the expressions relating various universal quantities to each 
other. If, under = 2, we apply the transformation (13) to the magnetization components in 
Eq. (1), the free-energy expansion remains the same, provided M4, are replaced according 
to Eq. (14) while the sextic coupling constants are transformed as prescribed below: 
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Uq + — + — ^ Uq + qQ + VQ, 
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It means that the following relations between the effective coupling constants should hold: 



, , q6iu,v) V6(u,v) / 3 , , 3 , , 3 , 
Mu, V) H \ — = Uq{u + -V, -V) + qeiu + -v, -v) + Ve{u + -v, -f ), 



2 3 3 
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As one can check, the expansions (10), (11), and (12) do really satisfy these relations. 
Moreover, consideration of all the remaining limiting cases (see, e. g. Refs. [^,0) shows 
that these RG series are in accord with their counterparts obtained earlier for the Ising 
26| , |27[| and 0(n)-symmetric models. 



Let us use the expansions just found for estimation of the universal critical values of 
the sextic effective coupling constants for the cubic ferromagnets, i. e., in the physical case 
n = 3. To obtain the numbers of interest from the divergent (asymptotic) RG series, a proper 
resummation procedure should be applied. Here we use the Pade-Borel-Leroy resummation 
technique which demonstrates high numerical effectiveness both for simple [Ising and 0(n)- 
symmetric] models p3| , p^ , p8| , p9| and for complicated anisotropic systems preserving their 



internal symmetries (see, e. g. Ref. for detail). Since the expansions of quantities 



depending on two variables u and v are dealt with, the Borel-Leroy transformation is taken 
in a generalized form: 



oo 



/(«, v)=Y: %i^v^ = I e-H'Fiut, vt)dt, F{x, y) = ^ T^^^^^- (17) 



U 

To perform an analytical continuation, the resolvent series 



i+] + h)\ 



oo n n-l „L n-l 

F(x,y,A) = EA"EE , (18) 

n=0 /=0 m=0 "'■ 

is constructed which is a series in powers of A with the coefficients being uniform polynomials 
in M, V and then Pade approximants [L/M] in A at A = 1 are used. 

With the four- loop RG expansions in hand, we can construct, in principle, three differ- 
ent Fade approximants: [2/1], [1/2], and [0/3]. To obtain proper approximation schemes, 
however, only diagonal [L/L] and near-diagonal Pade approximants should be employed 
pl| . That is why further, when estimating Wg, gg, and Wg, we limit ourselves to the ap- 
proximants [2/1] and [1/2]. Moreover, the diagonal Fade approximant [1/1] will be also 
dealt with although this corresponds, in fact, to the use of the lower-order, three- loop RG 
approximation. 

The algorithm for estimating of the universal critical values of sextic effective couplings is 
as follows. Since, in fact, we have the power series for the ratios Ru = Uq/u^, Rg = q^/iuv), 
and Rv = vq/v^ rather than for uq, qq, and vq themselves, we work with the RG series 
for these ratios. They are resummed in three different ways based on the Borel-Leroy 
transformation and the Pade approximants just mentioned. The Borel-Leroy integral is 
evaluated as a function of the shift parameter b under u = u*, v = v*. For the cubic 
fixed-point coordinates the values given by the six-loop RG analysis at n = 3 are adopted 
p!3| : u* = 1.321, V* = 0.117. The optimal value of b providing the fastest convergence of 
the iteration scheme is then determined. It is deduced from the condition that the Fade 
approximants employed should give, for b = bopt, the values of Ru {Rq, Rv) which are as 
close as possible to each other. Finally, the average over three estimates for R^ {Rg, Rv) is 
found and claimed to be a numerical value of this universal ratio. 

To demonstrate how such a procedure works, in Table 1 the results of the corresponding 
calculations are presented. The empty cells in this table reflect the fact that for some values 
of the shift parameter b Fade approximant [1/2] turns out to be spoiled by the "dangerous" 
poles, i. e., by the poles at positive or small negative t. As one can see, for Uq (an integer) bopt, 
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providing a maximal closeness of the estimates given by three working Pade approximants, 
is equal to 2, while for and fg, bopt = 3. So the results of our four-loop RG analysis are 
as follows: 

ul = 0.842, q; = 0.175, = 0.0108. (19) 

How close to the exact universal values of sextic couplings may these numbers be? Ob- 
viously, the accuracy of our estimates is limited both by the accuracy achieved in the course 
of evaluating the cubic fixed-point coordinates u*, v* themselves and by the speed of con- 
vergence of the iteration procedure employed. The influence of the latter factor may be 
characterized by the sensitivity of the numerical results with respect to variation of b. As 
is seen from Table 1, the values of Ug, q^, and fg, averaged over the working Pade approxi- 
mants, vary by no more than ±0.005, ±0.002, and ±0.0001, respectively, when b runs from 
to 20. On the other hand, these error bars are too small to be referred to as realistic ones. 
Their conservative counterparts may be deduced, accepting that the true universal values of 
renormalized sextic couplings should lie within the intervals of the variation of the estimates 
given by the leading Pade approximant, namely, by the higher-order pole-free approximant 
[2/1]. These intervals are easily extracted from Table 1 and are as follows: ±0.03 {uq), 
±0.01 (gg), ±0.001 (fg). Contrary to the aforementioned tiny ranges, the error bars thus 
determined look quite realistic. Moreover, they are found to be large enough to cover the 
inaccuracy produced by the uncertainty of the cubic fixed-point location. Hence, we adopt 
these error bars as the final ones. 

With the universal critical values of the quartic and sextic couplings in hand, we are able 
to estimate the reduced anisotropy of the sixth-order nonlinear susceptibility. Combining 
Eqs. (4), (5), and (19), we obtain 

5^^) = 0.102 ±0.02. (20) 



IV. EIGHTH-ORDER NONLINEAR SUSCEPTIBILITY AND DISCUSSION 

Let us discuss the role played by the sextic coupling constants in forming the magnitudes 
of xl^l Xd \ and 5(6). It may be shown that, in fact, their contrubutions to the sixth- 
order susceptibility and the anisotropy parameter are very small. Indeed, neglecting uq, 
qe, and vq in Eq. (4) changes the critical value of ^^^^ by about 4%. In other words, 
the anisotropy of x^^^ near Tc is fixed practically by the quartic coupling constants only. 
A similar situation should take place for the higher-order (eighth-order, tenth-order, etc.) 
susceptibilities. To demonstrate this, we consider the eighth-order susceptibility x^^^ fo^' the 
isotropic (Heisenberg) model. As x*^^'* and x^^\ it may expressed via the effective coupling 
constants 

= -40320^ (96ul - 24:U4Ue + us] . (21) 

177,9 V / 

Up to now, the universal critical value of the octic coupling constant was determined using 
the field-theoretical RG approach in three dimensions [^], the biased e expansion technique 



3^ , and the exact renormalization group machinery ||3^. The methods employed lead to 



the estimates = 0.168, Wg = 0.36, and = 0.145, respectively. Although these numbers 
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are considerably scattered, all of them are several times smaller than ul = 0.794 (Ref. p2[ ) 
and Mg = 0.951 (Ref. |^). It means that, because of the big numerical coefficients of 
the first two terms in brackets of Eq. (21), the contribution of ug to x^^"* is negligible at 
criticality. Since the cubic fixed point is located very near the Heisenberg one at the RG flow 
diagram, the same conclusion is valid for the 3D cubic model. This enables us to estimate 
the critical anisotropy of x^^^ without the calculation of the universal values of the octic 
coupling constants. Making in Eq. (21) the substitutions M4 ^ M4 + V4, Uq ^ Uq + Qq + Vq 
and U4 ^ U4 + ^4/3, uq + qe/3 + vq/9 with ug omitted, we obtain expressions for x^c^ 

(8) 

and Xd 5 respectively. Using then the known coordinates of the cubic fixed point and the 
critical values of the sextic coupling constants (19), we obtain the reduced anisotropy of the 
eighth-order susceptibility: 



.(8) _ v(,^)| 



^(8) ^ \Xc ^^Xd I ^ Q_^44 ^ g_g4_ (^22) 



The estimates (6), (20), (22) for S^^\ 5^^\ and S^^^ at criticality do not look too small to 
prevent the detection of the universal cubic anisotropy in physical and computer experiments 
or by means of a thorough analysis of the high-temperature expansions P^ , p5 |. It is well 



known, however, that the approach of the universal critical regime is controlled by the tiny 



exponent u ~ 0.01 |13], making both an experimental study and simulation of the true mode 



of critical behavior rather difficult |3^. We believe that the results obtained in this paper 



will help those searching for the cubic asymptotic regimes in real systems and lattice models 
to choose a proper strategy and the quantities to be measured. 



V. CONCLUSION 

To summarize, we have evaluated the nonlinear susceptibilities of the forth, sixth, and 
eighth orders for the 3D cubic ferromagnet and estimated the parameters S^'^\ S^^\ and 
5^^^ characterizing their anisotropy at the critical point. In the course of this study, we 
have calculated the four- loop RG expansions for the effective coupling constants Uq, q^, and 
Vq entering the small-magnetization free-energy expansion and found, resumming the RG 
series by the Pade-Borel-Leroy technique, the universal critical values of the renormalized 
sextic couplings at the cubic fixed point. The values of the anisotropy parameters (5*^*^ have 
turned out to be appreciable at criticality, indicating that the cubic (non- Heisenberg) critical 
behavior should be, in principle, detectable in physical and computer experiments. 
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TABLES 



TABLE I. Numerical estimates for the universal values of sextic effective coupling constants 

uq, qq, and vq obtained from the four-loop RG expansions (10), (11), and (12) resummed by the 
Pade-Borel-Leroy technique using approximants [2/1], [1/2], and [1/1]. The empty cells are due to 
the "dangerous" poles spoiling the corresponding Fade approximants. 



b 





1 


2 


3 


5 


10 


20 




[2/1] 


0.8654 


0.8514 


0.8427 


0.8368 


0.8294 


0.8206 


0.8144 




[1/2] 


0.8345 




0.8430 


0.8384 


0.8356 


0.8328 


0.8309 




[1/1] 


0.8097 


0.8295 


0.8400 


0.8466 


0.8543 


0.8625 


0.8679 


Q6 


[2/1] 


0.1834 


0.1791 


0.1765 


0.1747 


0.1725 


0.1699 


0.1680 




[1/2] 


0.1729 


0.1714 




0.1738 


0.1737 


0.1729 


0.1723 




[1/1] 


0.1653 


0.1707 


0.1736 


0.1754 


0.1775 


0.1797 


0.1812 


^6 


[2/1] 


0.01168 


0.01128 


0.01103 


0.01086 


0.01065 


0.01041 


0.01024 




[1/2] 


0.01062 


0.01059 




0.01090 


0.01072 


0.01064 


0.01060 




[1/1] 


0.00989 


0.01035 


0.01060 


0.01075 


0.01093 


0.01111 


0.01124 
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